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Abstract 

The ASCA resuhs of the starburst galaxy M82 are presented. The X-rays in the 0.5-10 keV band 
exhibit a thin thermal spectrum with emission lines from highly ionized magnesium, silicon, and sulfur, as 
well as a hard tail extending to higher than lOkeV energy. The soft X-rays are spatially extended, while 
Oh' the hard X-rays show an unresolved point-like structure with possible a long-term flux variability. The 

o 

5-H ' flux ratio of the emission lines and the spatially extended structure in the low-energy band indicate that at 

^ ' least two-temperature thin thermal plasmas are present. The abundances of the oxygen, neon, magnesium, 

(s^ silicon, sulfur, and iron in the thin thermal plasmas are found to be significantly lower than the cosmic 

C^ ' value. Neither type-la nor type-II supernova explosions can reproduce the observed abundance ratio. The 

origin of the unresolved hard X-rays is uncertain, but is probably an obscured low-luminosity AGN. 

Key words: Galaxies: individual(M82) — Galaxies: abundance — Galaxies: active — Galaxies: 
X-rays — X-rays: spectra 

1. Introduction 

Starburst galaxies comprise a sub-class of young galaxies, which exhibit violent star formation activities, mainly 
near to the center of the galaxy. The starburst activity may produce hot plasmas, due either to frequent supernova 
explosions or to strong stellar winds from young massive stars. The hot plasmas can partly escape to intergalactic 
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space as the galactic wind. Thus, the starburst galaxies would be good laboratories to study a possible interplay 

between star formation, heating of interstellar and intergalactic gas and chemical evolution over a wide range of 

the astronomical scales; from supernova and star clusters to the inter-galactic scale. X-ray observations should 

provide a key to study these issues, because the hot plasmas preferentially emit their energy in the X-ray band, 

giving information about the plasma temperature and density. Also, X-ray spectroscopy provides a direct method 

to determine the chemical composition. 

M82 has been known to be the most active starburst galaxy near the galaxy, and, together with its proximity 

of 3.25 Mpc distance, is placed to be the best-studied object (e.g., Rieke 1980). An X-ray image of M82 was 

firstly provided with Einstein HRI observation (Watson et al. 1984). They found, in addition to several unresolved 

sources, diffuse X-rays extending along the minor axis of the M82 galaxy. From the IPC data, Fabbiano (1988) 

determined the temperature of the extended emission to be 1.2-2.7 keV, while the MPC data revealed the existence 

of a high-temperature component, possibly associated with or near to the galactic nuclei. Although these authors 

failed to distinguish whether the diffuse emission is thermal or non-thermal, they interpreted that the X-rays are 

due to a galactic wind from the center of M82 at a rate of ~ 0.7 Mq yr^^. More detailed studies were made with 

ROSAT PSPC and HRI observations (Watson et al. 1994; Bregman et al. 1995). ROSAT HRI detected a point 

source at the galactic nucleus with Lx= 10**^ erg s^^. It even exhibited a flux variability on a time scale of about 

two weeks (CoUura et al. 1994). Strickland et al. (1997) extensively studied with ROSAT PSPC and HRI, and 

found that the diffuse emission is really thin thermal emission with a temperature of k^T— 0.6 keV at the center, 

and decreasing to kQT= 0.4 keV along the minor axis with absorption of N'ii=l x 10^^ cm^^. They also reported 

that the nuclear region exhibits two temperature components of k-QT= 0.76 keV with A^h^ 0.9 x 10^^ cm^^ and 

k-QTy 6.2 keV with Ah= 6 x 10^^ cm^^; the absorption of the hotter component is larger than those of the diffuse 

emission. Ginga found that the X-ray spectrum in the 2-10 keV band could not be fit with a power-law model, but 

could be nicely fit with thermal bremsstrahlung at a temperature of k^T — 5.75_Q'5ykeV absorbed by a cold column 

of A^H^ 8.9^3Q x 10^^ cm^^ (Tsuru 1992). The author also reported a marginal detection of the iron K emission 

line. However, the center energy was not determined well; we were thus unable to distinguish whether the hard 

X-rays are of thermal or non-thermal origin. All of these previous observations have been limited in either imaging 

capability (for the hard X-rays) or energy resolution. The first experiment with moderately high-energy resolution 

in the hard X-ray band was made with BBXRT. Petre (1993) reported, although only preliminary, the first detection 
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of K-shell transition lines of highly ionized magnesium, silicon, but no iron. However, the limited statistics prevent 

the author to obtain any detailed information. ASCA is the first X-ray satellite having wide energy-band imaging 

and spectroscopy capability with reasonably high statistics, and hence may be a suitable instrument for studying 

the nature of X-rays from M82. Since the image quality of the low-energy band is worse than those of the previous 

satellites, we put our emphasize more on the wide-band spectra than on the imaging aspects. Some preliminary 

reports concerning the ASCA observation have already been made by Tsuru et al. (1994), Awaki et al. (1996), Tsuru 

et al. (1996), and Ptak et al. (1997). Ptak etal (1997) focused on the SIS data and Moran, Lehnert (1997) reported 

results mainly from GIS and PSPC. This paper reports on the results of a combined study using both of GIS and 

SIS. 

2. Observation and Data Reduction 

M82 was observed with ASCA from 1993/04/19 22:22 to 04/20 07:58 using two sohd-state imaging spectrometers 
(SIS and SIS 1) with the 4CCD Faint and Bright modes and two gas-imaging spectrometers (GIS 2 and GIS 3) with 
the normal PH mode; both are separately placed on each focal plane of four identical thin- foil X-ray mirrors (XRT). 
Details concerning the instruments are described in Burke et al. (1991), Ohashi et al. (1996), Makishima et al. 
(1996), and Serlemitsos et al. (1995), while a general description of ASCA can be found in Tanaka et al. (1994). The 
data were screened with the standard selection criteria with XSELECT version 1.3 so as to exclude data affected 
by the South Atlantic Anomaly, earth occultation, and regions of low geomagnetic rigidity. We eliminated the 
contamination by the bright earth, removed hot and flickering pixels from the SIS data and rejected particle events 
from the GIS data using the rise-time discrimination technique. 

After these screenings, we obtained effective accumulation times of 18.4 and 16.0 ks for SIS and GIS, respectively, 
out of the total exposure of 22.5 ks. The counting rates of SIS and GIS were 1.3 counts s^^ SIS^^ and 0.86 
counts s^^ GIS^^, respectively. 

3. Analysis and Results 
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3.1. Spectral Analysis 

We collected X-ray photons from the 6' radius circle around the center of M82 image. The background spectrum 
of the SIS data was taken from a square region of 17'. 14 x 17'. 14, excluding the region r > 7'. 65 from the center 
of M82, while that of GIS was collected from the 7'. 65 < r < 18'. annulus around M82. After subtracting the 
background, we made a combined SIS spectrum (SIS and SIS 1 were added, SIS 0+1 hereafter). The data from 
GIS 2 and GIS 3 were also combined (GIS 2+3 hereafter). Figure 1 shows the thus obtained SIS 0+1 and GIS 2+3 
spectra. We used the XSPEC version 8.50 or 9.00 for the following spectral fits. 

The number of atoms per hydrogen for the cosmic metal abundance adopted in this paper are 0.0977, 3.63 x 10""*, 
1.12 X 10-'', 8.51 X 10-'', 1.23 x 10"'*, 3.80 x 10-^ 3.55 x 10-^ 1.62 x 10-^ 3.63 x 10~^ 2.29 x 10"^ 4.68 x 10~^ 
and 1.78 x 10"^ for He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni, respectively (Anders, Grevesse 1989). 

3.1.1. Low energy emission lines 

From figure 1, we can see several resolved emission lines above 1.2 keV and unresolved complex lines near and 
below 1 keV. We have determined the center energies and fluxes of the resolved lines by fitting with a model of 
narrow Gaussian lines on a power-law continuum in a limited energy band. The line energies and fluxes are shown 
in table 1 along with the most probable atomic-line candidates. The presence of emission lines from highly ionized 
magnesium, silicon and sulfur indicates that the X-rays are, at least partly, attributable to a thin thermal plasma. 

The line fluxes of H- and He-like atoms give the population of the ionization states, or the ionization temperature 
of the atom. The ionization temperature is a function of the electron temperature k^T and an ionization parameter 
Uct, where n^ and t are, respectively, the electron density and elapse time of the plasma with temperature of k^T . 
We note that at log(nct) > 12.5, the ionization temperature becomes nearly equal to the electron temperature 
(collisional ionization equilibrium or CIE in short), and that many astronomical plasmas are far from CIE; they are 
still in the ionizing phase, called non-equilibrium ionization, or NEI in short. 

From the silicon and magnesium lines listed in table 1, we made a two-dimensional parameter map of n^t and the 
electron temperature k-gT (figure 2) using the code of ionization in a plasma developed by Masai (1984, 1994). 

No overlapping region for silicon and magnesium was found, indicating that the plasma was not at a single 
temperature. We thus proceed to the spectral analysis with multi-temperature components, either in CIE or NEI. 
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3.1.2. Spectral fitting 

The SIS 0+1 and GIS 2+3 spectra were fitted simultaneously by a 2-temperature Raymond Smith thin thermal 
plasma model of CIE (Raymond, Smith 1977) at the temperatures of 1.1 and 0.5 keV, which were estimated from 
the line flux ratios given in figure 2 [see Si and Mg temperatures at log(nct) > 12.5]. We found a large excess in the 
data above the 3 keV bands. This immediately indicates that the spectra require another hard component whose 
temperature is above 10 keV. We thus conducted a simultaneous fitting to the SIS 0+1 and GIS 2+3 spectra with 
a three-temperature model (3-Raymond Smith or 3-RS): the soft, medium and hard components, hereafter. 

As shown in the next section, the soft component has an apparent extension. On the other hand, the emitting 
region of medium and the hard components are compact at the center region of the soft component. We thus assume 
the following model spectrum; 

iVn (Whole) x [Thermal(Soft) 

+ iVn (Medium) x Thermal (Medium) 

+ iVH(Hard) x Thermal(Hard)]. (1) 

We fit the spectra with a 3-temperature model using independent absorption columns for each component. Since 
the light metals, such as silicon and sulfur, are expected to be fully ionized at the high-temperature (hard) component, 
the abundances of Si and S in the hard component are not well constrainted. We therefore fixed the abundance ratios 
in the hard component to be solar ones. For the other components, the low- and medium-temperature plasmas, 
we let the metal abundances be free. With this model, the metal abundances of the medium and hard components 
are reasonably constraint, but not for the low-temperature component. No additional absorption for the medium 
component to the whole absorption (also for the soft component) is required with an upper limit of 3.2 x lO^^cm^^. 
We therefore assume that the metal abundances between the soft and the medium components are the same with 
no additional absorption for the medium component. This model gave a reasonable x^/d.o.f. of 444.8/372. The 
best-fit results are given in figure 1, 3 and table 2. We note that the metal abundances of the medium components 
are very similar to those obtained using the previous model, in which the abundances for both the soft and medium 
components are free. The relative normalization of GIS 2+3 to SIS 0+1 is free, and turns out to be 1.09, which is 
within the error on the cross-calibration of GIS and SIS. In short, we conclude, from a series of spectral fittings, 
that the abundances of the medium and hard components are well constrainted, and those of the soft component 
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are consistent with being the same abundances of the medium. Additional absorption is only required for the hard 

component. 

An alternative model, in which the hard component is replaced by a power-law model, gave a good fit with a 
X^/d.o.f. of 454.2/373. In this case also, additional absorption is required only for the hard power-law component. 
Thus, no essential difference was found between the non-thermal (power-law) and thermal (RS) models for the 
hard component; hence, whether the hard component is thermal or non-thermal remains uncertain. The best-fit 
parameters are given in table 3. 

Moran, Lehnert (1997) also reported that the X-ray spectrum of M82 can be described by three components: two 
thin thermal plasmas and a power law dominating at the high-energy band. However, this result is different from 
ours; the large cold absorption column to the medium component, or the "warm thermal component" in Moran, 
Lehnert (1997), is not required in our fitting. This discrepancy would be attributable to the assumption on a 
metal abundance. Moran, Lehnert (1997) assumed solar abundances for the two thermal components. This model, 
however, does not fit the SIS line structures with higher energy resolution than those of GIS and PSPC. 

3.1.3. Upper limit on iron K emission line 

Since the iron K-shell emission line provides key information concerning a plasma diagnostics of a temperature 
higher than a few keV, we made a spectral fitting with a thermal bremsstrahlung model plus a Gaussian line around 
6-7 keV in a limited energy band of 4-8 keV. Since the energy band was selected to be reasonably narrow, the 
result did not depend on the spectral model of the continuum component. We found no significant emission line. 
The upper limit is shown in figure 4 as two-dimensional (Line energy-equivalent width) contours. Two-dimensional 
contours determined with Ginga are also given (Tsuru 1992). 

3.2. Imaging Analysis 

Figure 5 shows the SIS images in the 0.5-0.8 keV, 1.2-1.8 keV and 3-10 keV bands: the energy bands which 
represent the soft, medium and hard components, determined from the 3-RS fitting (see figure 1). 

We found that the 0.5-0.8 keV band X-rays show significant spatial extension, particularly along the minor axis 
of the optical image of M82. The radial profile in the range of r = 0'.5-5' can not be represented by the point spread 
function (PSF in short) with x^/d.o.f= 35.9/12. This, together with the spectral information, indicates that the 
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soft component is a diffuse thermal plasma extending larger than the optical image of M82. 

The radial profile of the 3-10 keV band, on the other hand, is consistent with the PSF (x^/d.o.f= 8.8/12). The 
upper limit on the spatial extension is 0'.5 in radius. Thus, the hard component is point-like within the ASCA 
resolution. 

The radial profile of the 1.2-1.8 keV band image is not fitted with the point spread function with x^/d.o.f= 27.6/12. 
Therefore, the medium component is extended, but its extension may be smaller than that of the soft component. 
The peak positions of the 1.2-1.8 keV and 3-10 keV bands agree with each other within the spatial resolution. The 
0.5-0.8 keV band image is more complex, as is shown in the closed view (figure 6). A couple of peaks around the 
galactic center are seen in the 0.5-0.8 keV image. However, neither peak agrees with the peak of the 1.2-1.8 keV 
and 3-10 keV bands. 

4. Discussion 

We found that the X-rays of M82 have a very complex structure, comprising extended soft and medium- 
temperature plasmas and point-like hard X-rays near to the center of M82. Based on the ASCA results, and 
referring to previous X-ray observations, we try to give possible implications on the nature of M82 X-rays. 

4-.1. Hard X-Ray s 

4-. 1.1. Time variability 

The spatial extension of the hard X-rays near to the center of M82 is constrainted to be less than 0.5 kpc (0'.5) 
radius. Moran, Lehnert (1997) already reported that no short time variability was detected in the light curve of 
GIS 2+3 with an energy band of 3.5-10 keV. We also detected no short time variability in SIS 0+1. 

The hard component should be the same as those observed with previous non-imaging hard X-ray instruments, 
from the first detection of Uhuru to recent Ginga scanning observations. Since most of these previous observation 
might have suffered from possible contamination of M81, lying about 1 degree away from M82, we compare them 
with the results from relatively small field-of-view experiments of EXOSAT ME and scanning observation of Ginga 
(Schaaf et al. 1989; Tsuru 1992). Both of them safely avoided M81 contamination. When calculating the X-ray flux 
for each experiment, we adopted the best-fit spectral model to each observational spectrum. The results concerning 
the X-ray fluxes in the 2-10 keV band are, (2.3 ± 0.2) x 10"i\ (3.4 ± 0.1) x 10"i\ 1.9 x 10"" erg cm'^ s^^ 
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for EXOSAT ME, Ginga and ASCA, respectively. Assuming that M81 contributed 15% of the flux obtained from 

EXOSAT ME observation, the M82 flux observed with EXOSAT ME becomes (2.0 ± 0.2) x 10"" erg cm"^ g-i 

(Schaaf et al. 1989). 

We see a significant variation of about a factor of 2 in the X-ray luminosities. Since 90% of the flux in the 2-10 kcV 
band derived from ASCA was emitted from the hard component, we conclude that the hard component of M82 has 
been variable with a long time scale of years. Adding to this, a significant spectral change from Ginga to ASCA 
is seen. The X-ray spectrum at the time of the Ginga observation was not fit with a power-law model, but was 
nicely fit with thermal bremsstrahlung at a temperature of k^T — 5.75_Q'5y keV, absorbed by a cold column of 
A^H= 8.9110 ^ 10^^ cm"2 (Tsuru 1992), which is different from the spectral shape obtained with ASCA. 

These two results concerning the flux and spectral shape indicate that the source of the hard component is a time 
variable. The time scale gives the upper limit on the size of the hard X-ray emitter to be less than 1 pc, which is a 
far more severe constraint than that obtained from imaging analysis. 

The hard component might be a summation of the unresolved galactic point sources. However, at least half 
of the flux in the 2-10 keV band observed with Ginga was time variable with a timescale of several years, which 
immediately indicates that at least half the time variable was due to a point source (Tsuru 1992). 

4- 1-2. Comparison with ROSAT observation 

The spectrum of the nuclear region obtained with ROSAT PSPC is reproduced by a two-thermal bremsstrahlung 
model of k-QT> 6.2 keV and k-QT= 0.76 keV (Strickland ct al. 1997). This may be consistent with our interpretation: 
a hard X-ray source at the M82 center and diffuse soft and medium X-rays around the center. 

Near to the M82 center, ROSAT HRI found three point sources; among them, source 8 is the brightest, by a factor 
of 25, than any other point sources (Strickland et al. 1997). Therefore, this HRI source could be a possible candidate 
of the ASCA hard X-ray source. The position error of the ASCA image, however, is as large as about 0'.5 and no 
reference X-ray star in the ASCA field-of-view is available. Furthermore, as we already noted, the peak of the hard 
and medium components does not agree with that of the soft component. Thus, identification of the hard X-ray 
source to source 8 is debatable. Nevertheless, we discuss the possibility of the hard component to be source 8. From 
the best-fit spectral model parameters given in table 2 and HRI spectral response, the flux of the hard component 
obtained from our observation was converted to a HRI counting rate of 0.018 counts s"^. This value is 2-3 times 
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smaller than the observed ROSAT HRI flux of source 8, 0.050 counts s^^, but within the variable range obtained 

during the HRI observation (CoUura et al. 1994). 

4-.1.S. Origin of the hard X-rays 

Rieke et al. (1980), Schaaf et al. (1989), and Moran, Lehnert (1997) argue that an inverse-Compton X-ray emission 
scattered from infrared photons is an origin of the hard X-rays. However, the size of the emitting region, evaluated 
from the observed long-term variability, is estimated to be only ~ 1 pc, far smaller than the previous estimation by 
an order of two or three. Thus, the contribution from this process to the total X-ray emission could be ignored. 

The brightest ROSAT HRI source 8, can be identified with a bright radio source, either 41.5-1-597 or 41.95+575 
(Strickland et al. 1997). The radio source 41.95+575 may be attractive as a counterpart of the time variable ASCA 
hard X-ray source or ROSAT HRI source 8, because its radio flux is variable with a decreasing rate of 9% per a 
year. Adding to this, its synchrotron self-absorption frequency is also decreasing and the source shows a shell-like 
structure with a diameter of 0.3 pc (Muxlow et al. 1994). All of these results led Kronberg et al. (1985) to propose 
that 41.95+575 is a type-II supernova, exploded in a very dense region with no ~ 10^ cm^^. Terlevich (1994) also 
modeled that 41.95+575 is an SNR of the type-II supernova which occurred in ^ 1955. This putative young SNR 
(41.95+575) may produce hard X-rays, as is found in other young supernova remnants, such as SN 1993J in M81 
and SN 1978K in NGC 1313 (Kohmura et al. 1994; Petre et al. 1994). The initial X-ray luminosity of theses SNRs 
are about 10^° erg s^^ (SN1993J), which is more or less comparable to the hard X-ray flux of M82. However, unlike 
the above sample of young SNRs, the hard X-ray luminosity of M82 has decreased from Ginga to ASCA by factor 
of 2, while the temperature has increased from 5.75 keV to 14 keV. Furthermore, no strong iron-K line was found. 
Thus, the X-rays from a possible young SNR (the radio source 41.95+575) is unlikely for the origin of the hard 
component. 

One may argue that the hard source would be a super-Eddington source, which is often found in some spiral 
galaxies. However, the luminosity of the hard component of 4 x 10^° erg s^^ is far larger than that of any super- 
Eddington source (10"^^ to 10^° erg s^^, see e.g. Okada et al. 1994). 

Another possible counterpart would be 41.5+597, which is highly variable in the radio band; it was ^ 7 mJy in 
1981 and decreased by ~ 100% within a year (Kronberg, Sramek 1992). An HST observation showed that there is a 
very bright spot at the position 41.5+597 (region E in O'Connell et al. 1995). In this case, the most likely scenario 
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for the hard X-rays is a low-luminosity AGN (LLAGN hereafter) behind the dense gas near to the center. The 

absorption of the hard component is A^h=(2-3) x 10^^ cm^^. The A^h value estimated from the CO observations 

at the position of the hard component (i.e. at the source 8) is A'h=(4 6) x 10^^ cm^^ (Lo ct al. 1987; Nakai et al. 

1987). Assuming that about half of the A^h is located in front of the M82 center, the ASCA result is consistent with 

the scenario of the absorption due solely to the molecular clouds. 

Several ASCA observations show that some spiral galaxies are hosts of LLAGN, whose luminosities are 10^" to 

10"*^ erg s^^ (Ishisaki et al. 1996). ROSAT HRI detection of a short-time variability may also support this scenario 

(CoUura et al. 1994). One may argue that no clear AGN activity has been reported in M82. However, the large 

absorption column and the fairly complicated starburst activities would make it hard to detect any LLAGN activity 

at other wavelengths. 

4-. 2. The Soft and Medium Temperature Plasmas 

We found two-temperature extended plasmas from M82. Since the soft component is more extended than the 
medium component, the absorbing column for these components is larger than the galactic value of 4.3 x lO^^cm^^, 
it needs extra absorption near to M82 (Stark et al. 1992). The metal abundances of these hot plasmas were also 
constraint to some extent. We will discuss these issue separately. 

4-2.1. The X-ray images and relevant absorption 

We see that the counting rate in the 0.5-0.8 keV band at the peak of the other two band images in figure 6 is low. 
This would be due to the absorption of cool matter around the M82 center. Comparing the 0.5-0.8 keV band image 
with the CO and H I images, we found that the low-brightness region in the 0.5-0.8 keV band corresponds to the CO 
and H I disk (Lo et al. 1987; Nakai et al. 1987; Weliachew et al. 1984). The A'h value at the disk estimated from the 
CO and H I is (4-6) x 10^^ cm^^, which is large enough to absorb soft X-rays. Therefore, a simple idea is that the 
diffuse soft X-ray emitting plasma has a flux peak at the center as well as a component of the medium-temperature 
plasma, but highly absorbed by molecular clouds around the center. 

The soft and medium components extend more than the thickness of the molecular disk of ~ ±100 pc. Its average 
column density between \z\ — 100-500 pc is A'h^ 6 x 10^^ cm^^ (Nakai et al. 1987). Then, for the relevant X-rays, 
the mean absorption column by the molecular clouds would be A^h^^ 3 x 10^^ cm^^, half of the total A^h. The ASCA 
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observation shows the absorption due to the M82 cool gas to be A^h^ (2.6^q'4) x 10^^ cm^^ after subtracting the 

galactic absorption from the total absorption of (S.O^q^) x 10^^ cm^^ , which agrees well with the results of CO 

observations. 

4-2.2. Metal abundances 

It is expected that a large number of sequential supernovae, mostly type-IIs, explode as the result of starburst 
activity, synthesize a huge amount of metals and cause strong emission lines in X-rays. However, the observed 
metal abundances are all below the cosmic values, especially the abundance of oxygen, the most important product 
in type-II supernovae, is found to be extremely low. Suppose that all the energy of the hot plasma (i.e. the soft 
and medium components) is supplied by supernova explosions, and the origin of the plasma gas is a mixture of the 
supernova ejecta and interstellar gas; then, the thermal energy of the hot plasma per one supernova explosion can 
be written as 

3M^+M^^^^ < Es., (2) 

2 /i?TT,H 

where Mcj and Mam are the masses of the ejecta and the mean mass of the ambient gas heated by a single supernova. 
The ratio of the mass of each metal (M*) to the total mass ( Mtotai) is expressed as 

M* _ M* + M:^ > M* 



Mtotai Mej + Mam '^ Mej + Mam ' 

where Mam and M*- are the mass of each metal in the ambient, and are synthesized by the supernova explosion. 
From equations (2) and (3), we obtain 



Mtotai ^ 2 ""J limn Esn 



4.9 X 10-3 ( 4^^' ' 
lO^^erg 



kT,T \ / M 



cj 



-^ 



(4) 



1 keVy VI MQj 

We adopted the results of m„ = 50 Mq by Tsujimoto et al. (1995) as M*:, in which the synthesized masses are 1.8, 
0.23, 0.12, 0.12, 0.041, 0.0080, and 0.091 Mq for O, Ne, Mg, Si, S, Ar, and Fe, respectively. Using these values, 
equation 4 and the cosmic abundances in Anders and Grevesse (1989), the lower limits on the metal abundances are 
estimated to be 0.71, 0.51, 0.72, 0.66, 0.43, 0.31, and 0.18 for O, Ne, Mg, Si, S, Ar, and Fe, respectively, based on 
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the assumption of _Esn = 10"''^ erg and k^T— 1 keV. The observed abundances of sulfur and argon agree with this 

assumption, but those of oxygen and iron are smaller than the results based on it. 

One may argue that the low abundances can be explained by the remaining dust in the hot gas. The number of 

oxygens synthesized by type-II supernovae is about 10 times larger than those of magnesium, silicon, sulfur, and iron. 

The number density of oxygens observed in M82 is still higher than the sum of the number densities of magnesium, 

silicon, sulfur, and iron by factor of 1.6. Even in the extreme case that oxygen, magnesium, silicon, sulfur, and iron 

would be confined in the dust as an oxide, a large fraction of the oxygen could be free from the dust. Neon is also a 

free atom. Thus, the dust-condensation scenario can not explain the low metalicity of the oxygen and neon in M82. 

4-2.3. Relative ratios of the metal abundances 

Reserving the potential problem of the absolute metal abundances, we further discuss the relative abundance. We 
compare the metal abundance ratios of M82 and those of type-II supernova in figure 7. The metal abundances of 
type-II supernova are integrated and averaged over the possible SN mass range with a reasonable mass function 
(Tsujimoto et al. 1995). For a comparison, we also plot the values of type-la supernovae in figure 7, which are 
adopted from model W7 in Nomoto et al. (1984) and Thielemann ct al. (1986). From figure 3, we can see that 
the abundance of the M82 gas peaks at the medium elements, such as silicon and sulfur, which are found neither 
in type-II nor typc-Ia supernova. This does not agree with the scheme that metal synthesise due to supernova 
explosions contributes to the metal abundances in M82. 

Since the values of the abundance ratios of light elements in M82 are between those of type-la and type-II 
supernova, one may argue that they can be explained by a mixture of the two types of supernova, although the 
type-la supernovae may be minor compared to type-II in a starburst galaxy. However, the abundance ratio of iron 
to silicon in M82 is well below both that of type-la and type-II supernova. Accordingly, we have no solution which 
satisfies both the abundances of iron and medium elements with this scenario. 

We note that although other starburst galaxies show similar low metal abundances to that of M82 (e.g. ASCA 
observation of NGC 253; Awaki et al. 1996), starbursts associated with the Seyfert activity (e.g. NGC 1068; Ueno et 
al. 1994) exhibit normal metal abundances. This may imply that starburst activities, or relevant metal abundance, 
may differ from object to object. 

One of candidates for the origin of metals in the intra-cluster medium (ICM hereafter) is ejecta from starburst 
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activities during the early phase of galaxy evolution (e.g. Sarazin 1988 and references therein). Although M82 is 

a nearby galaxy and is not a proto-galaxy, it exhibits typical starburst activity and is a unique target to address 

this issue at this moment. The ASCA observations of four bright clusters of galaxies reveal that the average metal 

abundances of oxygen, neon, silicone, sulfur, and iron are 0.48, 0.62, 0.65, 0.25, and 0.32, respectively (Mushotzky 

et al. 1996). The observation of M82 shows lower abundances for all of the metals. In particular, oxygen and iron 

abundances of M82 are lower than those of ICM by a factor of 8-10. This suggests that the ejecta from starburst 

activity would not be the origin of ICM. 

From the large number of samples obtained with ASCA, Matsumoto et al. (1997) showed that the metal abundances 

of the hot inter stellar medium (ISM hereafter) in elliptical galaxies are about 0.3. The average metal abundances 

of hot ISM in three bright elliptical galaxies obtained with ASCA are 0.37, 0.33, 0.41, and 0.40 for oxygen, silicon, 

sulfur, and iron, respectively (Awaki et al. 1994). The abundances of silicon and sulfur in M82 are comparable to 

those of the elliptical galaxies, but those of oxygen and iron are less. 

4-. 3. Non Equilibrium Ionization 

Since the NEI model predicts a higher electron temperature than the line-emitting ionization temperature, it may 
be conceivable that low-temperature NEI model can reproduce the overall spectra. We thus tired to simultaneously 
fit the SIS 0+1 and CIS 2+3 spectra with a 2-temperature ionization non-equilibrium plasma model (2-NEI model 
in short). Since no NEI code based on the Raymond Smith thermal plasma model is available in the public XSPEC, 
we installed a NEI model based on the code of Masai (1984, 1994). We note that the total number of free parameters 
between 3-RS and 2-NEI are nearly the same as that of each other. This 2-NEI model gives a nice fit with x'^/d.o.f. 
— 441.3/390, almost equal to the case of 3-RS. The best-fit results are given in table 4. 

The small size estimated from the observed time variability of a couple of years gives a severe constraint on the 
NEI model. The emission measure n^T^ is 2 x 10^^ cm~'^, where Uc and V are the density and volume of the plasma. 
If the size of the emitting region is as small as 1 pc, then n^ is larger than 5 x 10'^ cm^'^. Therefore, the observed 
Uct value of 9 X 10^ s cm^^ constrains the age of the plasma to be smaller than three weeks. This extremely short 
time scale led us to conclude that the NEI plasma is unrealistic. Thus, we safely conclude that the origin of the 
hard X-rays is not a NEI plasma, but either a CIE plasma or a non-thermal source. 



14 T.G.Tsum ct al. [Vol. 49, 

5. Conclusion 

The X-ray structure of M82 composes of a point-like, heavily absorbed hard X-ray source with a kBT=14 keV 
temperature or power-law with photon index of 1.7, and extended two-temperature plasmas with kBT=0.3 keV and 
0.95 keV. 

We found a long-term variability and spectral change in the hard X-rays, and inferred that the most likely origin 
is due to an obscured low-luminosity AGN. 

The two-temperature plasmas would be produced by the starburst activity of M82. However, the observed low 
metal abundances in M82 gas are not easily explained by any mixture of supernova explosions. 

A two-temperature NEI model can reproduce the observed spectrum, but is physically unrealistic. 

The authors would like to thank Drs. R. Petre, R. Mushotzky, T. J. Ponman, R. J. Terlevich, D. K. Strickland, 
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supplied his NEI code to us. We wish to thank the ASCA software team, who developed the analysis software which 
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Figure Captions 

Fig. 1 SIS 0+1 and GIS 2+3 spectra obtained within a 6'.0 radius from the center of M82 and the best-fit 3-RS model. 

Fig. 2 Two-dimensional parameter map of riet and feT in the case of the NEI model. 

Fig. 3 Metal abundances from the 3-RS model fit. 

Fig. 4 Confidence contour maps between the iron K-emission line center energy and its equivalent width obtained with 
Ginga and ASCA. 

Fig. 5 Three bands images of SIS smoothed with a = 0'.212. The contour levels are 0.00969, 0.0145, 0.0218, 0.0327, 
0.0491, 0.0736, 0.110, 0.166, 0.248, 0.373 counts arcmin"^ s"^ The dashed lines show gaps among the CCD chips in SIS 0. 

Fig. 6 Closed version of the 0.5-0.8 keV image of SIS smoothed with a = 0'.159. The contour levels are 0.0194, 0.0242, 
0.0291, 0.0339, 0.0388, 0.0436, 0.0485, 0.0522, 0.0630, 0.0727, 0.0921 counts arcmin"^ s"\ The dashed lines show the gaps 
among the CCD chips in SIS 0. The crossed mark indicates the peak position of the medium and hard components. 

Fig. 7 Relative metal abundaces synthesized by type-la and type-II supernovae. The data have been adopted from Tsuji- 
moto et al. (1995), Nomoto ct al. (1984), and Thielemann et al. (1986). 
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Table 1 The line energies and fluxes*. 


Line Energy 


photon flux 


ID 


(keV) 


(10"'"' photon s^i cni-2) 




1.351 ±0.011 


7.5 ±2.3 


Mg Ka (He-like) 


1.457 ±0.043 


1 O+2.0 
^••^-1.3 


Mg Ka (H-like) 


1.864 ±0.008 


o 1+2.3 
8.1_2.4 


Si Ka (He-like) 


2.006 ±0.012 


5.8_2^o 


Si Ka (H-like) 


2.457 ±0.016 


2 5+2-0 

^•^-1.9 


S Ka (He-like) 



* All the errors including upper limits are described at 90% confidence limits. 
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Table 2 The three Raymond- Smith model* 



Parameter 




Component 








Soft Medium 


Hard 


Lx MIO^" erg s" 


') 


1.6 1.4 


3.6 


A'Hwhoio (10^^ cm" 


-') 


3.olS:f 




iVHhard (lO^l Cm- 


-') 




18.9lti 


ksT (keV) 




0.32+_rol 0.95tS:^^ 


i3.9f^:^ 


Metahcity (Cosmic)'^ 


N 


< 0.15 









0.063t°:°?^ 






Ne 


U--LO-0.03 






Mg 


n oc+o.o? 


0.15^ 




Si 


0.40troI 






S 


0.47«;1? 






Ar 


< 0.51 






Fe 


o.049t°o;gj^ 




xVd.o.f. 


444.8/372 



* All the errors including upper limits are described at 90% confidence limits. 

a: Unabsorbcd luminosity in the 0.5-10 keV band. 

b: Galactic absorption of 4.3 x 10^°cm^^ is included (Stark et al. 1992). This parameter is common among the three 

components. 

c: Relative abundance among the elements are fixed to be cosmic values. 

d: Relative metal abundance ratios in the hard components is fixed to the cosmic values. 
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Tabic 3 The two ] 


laymond-Smith and power-law model*. 


Parameter 




Component 






Soft Medium 


Hard 


Lx ^ (10« erg s" 


') 


1.6 1.4 


4.2 


A^Hwhoic (lO^i cm 


-') 


o n+0.5b 




^Hhard (10^1 Cm- 


-') 




23.itl;g 


fcsT (keV) or I 




r)oi+0.04 A q^+005 


1.70 ±0.09 


Metalicity (Cosmic) 


N 


< 0.15 









n ofii+°o35 






Ne 


n. 1^ + 0.05 






Mg 


9fl+"-07 






Si 


0.40_oog 






S 


U-45_Q IQ 






Ar 


<0.48 






Fe 


o.048t°;°J^ 




xVd.o.f. 


454.2/373 



* All the errors including upper limits are described at 90% confidence limits. 

a: Unabsorbed luminosity in the 0.5-10 keV band. 

b: Galactic absorption of 4.3 x 10^°cm^^ is included (Stark et al. 1992). This parameter is common among the three 

components. 
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Table 4 The two non- 


equilibrium ionization model*. 


Parameter 




Component 
Soft Hard 


Lx '^ (10''° erg s" 


-') 


1.7 10.0 


iVn whole (10^1 cm 


-') 


1.9 ±0.2'^ 


^Hhard (10^1 cm- 


-') 


14.81^:? 


feT (keV) 




0.65t°;°^ 13.9lli 


log(nt) (s cm^"^ 


) 


> 11.6 9.96to:?{ 


Metalicity (Cosmic) 


N 


< 0.41 







n 29+°°^ 
u.zy_Q Qg 




Ne 


n 92+009 
"•^^-0.07 




Mg 


n 17+0-03 
"•^'-0.03 




Si 


n iq+0-04 




S 


A oc+0.05 
"•^0-0.06 




Ar 


<0.49 




Fe 


0.075t°;^^l 


xVd.o.f. 


441.3/390 



* All the errors including upper limits are described at 90% confidence limits. 

a: Unabsorbed luminosity in the 0.5-10 keV band. 

b: Galactic absorption of 4.3 x lO^Ocm^^ is included (Stark et al. 1992). This parameter is common between the 

two components. 
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Fugure 1: 



Figure 2: 



Figure 3: 
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Figure 4: 
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Figure 5: 
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Figure 6: 



Figure 7: 
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